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Abstract 

We study within the light-cone path integral approach the induced gluon emission from 
a fast quark passing through a finite-size QCD plasma. We show that the leading log 
approximation used in previous studies fails when the gluon formation length becomes of 
the order of the length of the medium traversed by the quark. Calculation of the energy 
loss beyond the leading log approximation gives the energy loss which grows logarithmi- 
cally with quark energy contrary to the energy independent prediction of the leading log 
approximation. 



In resent years much attention has been attracted to the problem of the induced gluon 
radiation from fast partons in a hot QCD medium (for a review, see Jl|). It is of great 
interest in connection with the current experiments at SPS, RHIC, and future experiments 
at LHC on A + A collisions since jet quenching due to the parton energy loss can be a 
good probe of formation of a hot quark-gluon plasma (QGP). 

Evaluation of the gluon emission from a fast parton in a medium requires the under- 
standing of the non-abelian analogue of the Landau-Pomeranchuk-Migdal (LPM) effect 
|2|, |J. There are two approaches to the LPM effect in QCD: the so-called BDMS approach 
Q (see also 0, ||) based on the Feynman diagrammatic formalism, and the light-cone path 
integral (LCPI) approach developed in our paper |J (see also [0, || ||, [L(J]). The BDMS 
approach neglects the mass effects, and applies for large suppression of the radiation rate 
as compared to the Bethe-Heitler one. The LCPI approach applies for arbitrary strength 
of suppression. For large suppression these approaches are equivalent j|, |], |IlJ]. The 
probability of gluon emission in the BDMS and LCPI approaches is expressed through 
the solution of a two-dimensional Schrodinger equation with an imaginary potential. This 
equation describes evolution of the color singlet qgq system in the medium. The poten- 
tial is proportional to the cross section for scattering of the qgq system on a medium 
constituent. For the QGP the constituents can be modeled as Debye-screened colored 
Coulomb scattering centers |E| . 



In H the quark energy loss, AE, has been evaluated analytically treating interaction 
of the qgq system with the Debye-screened centers in the Leading Log Approximation 
(LLA) which is equivalent to the harmonic oscillator approximation for the Hamiltonian 
of the qgq system. For a quark produced inside a finite-size QGP the BDMS prediction is 

AE BDMS = — — v, (1) 

4 Ag 

where L is the length of QGP traversed by the quark, /i is the Debye screening mass, X g 
is the mean free path of the gluon in QGP, Cf is the color Casimir for the quark, and the 
factor v grows smoothly with L, at L ^> X g v ~ log(L/A 9 ). 

The energy independent AE ([I]) differs from that obtained recently by Gyulassy, Levai, 
and Vitev |13| . Calculating the Feynman diagrams for the single scattering (the first order 
(N = 1) in opacity) they have obtained 

C F a s L 2 /! 2 E 

AE GLV = — — log - . 2 

4 X g 

Since the AEbdms should include the N — 1 contribution the contradiction between 
(H|) and (§) at E — > oo seems to be surprising^. By now there has not been given any 
explanation of this fact, except the argument of the authors of Ref. [I3| that it can be 



connected with the neglect of the finite kinematic bounds in the analysis |4| . However, it 
is clear that it cannot be important at E — > oo. 



1 Strictly speaking, the derivation of the BDMS formalism given in Ref. H is valid only when the 
number of rescatterings is large. However, since the formulas obtained are equivalent to those of the 
LCPI M approach which is free from this restriction, it is clear that the BDMS prediction should contain 
the N = 1 term. 
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In the present paper we resolve the above puzzle of the discrepancy between the 
BDMS and GLV predictions. We demonstrate that the absence of the logarithmic energy 
dependence in ([!]) is connected with the fact that the LLA fails when the gluon formation 
length becomes of the order of L. In this case the spectrum is dominated by the N = 
1 scattering which simply vanishes in the LLA. We show that if one uses the actual 
imaginary potential the energy loss grows logarithmically with quark energy. However, 
the denominator in the argument of the logarithm is not the Debye mass as it is in (|2|). 

We will work in the LCPI formalism 0. The probability distribution of the induced 
gluon emission from a quark produced at z = can be written as [10 



dP 

dx 



dz n(z) 



dcrf f H f (x,z) 



dx 



(3) 



where x is the gluon fractional momentum, n is the number density of the medium, and 



dx 



Re J dp^*(p,x)a 3 {p,x)^ m (p, 



x,z) 



(4) 



Here a 3 is the cross section for interaction of the qgq system with a scattering center. The 
relative transverse separations in the qgq system are p gq - = (1 — x)p, p qq - = —xp . ^(p, x) 
is the light-cone wave function for the q —> gq transition in vacuum, and \P m (p, x, z) is 
the quark light-cone wave function in the medium at the longitudinal coordinate z (we 
omit spin and color indices). The wave functions (modulo a color factor) read 



P(x) 



ty m (p,x,z) = P{x) 



d 



d 



IS, 



r ) j d£exp (-^-j Ko{p,t\p',0) 



v/ 







is 



d 



dp'x dp 



v/ 



p=o 
K(p,z\p ',*-£) 



p=o 



(5) 
(6) 



where P(x) = iJ a s /2x[s g (2 — x) + 2s q x]/2M(x), s q ^ g denote parton helicities, /C is the 
Green's function for the two-dimensional Hamiltonian 



H(z) 



2M(x 




.n(z)a 3 (p,x) 



(7) 



and 



fC (p 2 , z 2 \p 1 ,zi) 



M(x) 



2ni(z 2 - zi) 



exp 



iM{x){p 2 - Pl ) 
2{z 2 - z x ) 



is the Green's function for the Hamiltonian (0) with v(p, z) = 0, M{x) = Ex(l — x), and 

x) + m 2 x. The gluon mass m g plays the role of 



L 



f 



2Ex(l 



x)/e 2 with e i 



infrared cutoff removing the contribution from long wave gluons which cannot propagate 
in the QGP. It is natural to take m g ~ /i. However, for large suppression which occurs at 
E — » 00 the parton masses can simply be neglected. 



3 



The three-body cross section can be written as fll4 



C 1 
0-3 (P, x) = 2^rM(l - X )P) + ^(p) - jp^{xp)} , (9) 

where Ca = N c is the octet color Casimir, 02 (p) is the dipole cross section for scattering 
of a qq pair on a color center. For the parametrization 02 (p) = C*2(p)p 2 the factor C*2 is 

Here CV is the color Casimir of the scattering center. In the region /) < which 
dominates the spectrum for strong suppression fllCf ) takes the form 



, , C F C T 0?/K ( 1 \ 

Ca P ~ 9 log — • 11 

The LLA consists in replacing C 2 (p) by C 2 (p e //), where p e // is the typical value of 
p. This seems to be a reasonable procedure since C%{p) has only a slow logarithmic 
dependence on p. Then <7 3 (p, a;) = Cz{x)p 2 , where C^{x) = C2(p e ff)A(x) with A(x) = 
[1 + (1 — x) 2 — x 2 /N 2 ]Ca/2Cf, and the Hamiltonian (^) takes the oscillator form with 
the frequency Q(x) = \J —iC 3 (x)n/ M (x) . The value of p e ff is connected with the gluon 
formation length, //, by the Schrodinger diffusion relation p 2 ^ ~ lf/2M. If is simply the 
typical scale of £ in (5), (|) when the wave functions are substituted in ([|). 

Let us discuss the gluon emission at qualitative level. We begin by estimating p e jf 
and If. Let us first estimate these quantities for gluon emission from a quark in an 
infinite medium. We will denote them as p e ff and If. They should also be related by 
the Schrodinger diffusion relation. On the other hand, the absorption effects for the qgq 
system should become strong at the scale If. It means that IfnC 3 plff /2 ~ 1 . From these 



conditions one gets p e ff ~ [E a x(l — x)nC s ] 1 l i and l e ff ~ 2y E a x(l — x)/nC^. These 
estimates are valid when p e ff < 1/e and If <Lf. 

Now we turn to the gluon emission from a quark produced inside a finite-size medium. 
In this case in the high-energy limit qualitatively two different situations are possible. 
The first regime gets for the gluons with x such that If < L. In this case the finite-size 
effects play a marginal role, and p e ff ~ Peff- The spectrum can roughly be calculated 
using the effective Bethe-Heitler cross section for the infinite medium. We call this regime 
the infinite medium regime. The second regime occurs for the gluons for which If > L. In 



this case p e // ~ Pd(L), where pd(L) = \J L/2M is simply the diffusion radius on the scale 
of the quark path length inside the medium. In this regime the effective Bethe-Heitler 
cross section is chiefly controlled by the finite-size effects. We will call this regime the 
diffusion regime. Thus we can write for the above two regimes 

p eff ~ min(p eff ,p d (L),l/e). (12) 

Here we have taken into account that p e ff < 1/e. In terms of x the infinite medium 
regime occurs at x < 6 and (1 — x) % 5, and the diffusion regime gets at 5 < x < (1 — 5), 
where 

(13) 
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For the sake of definiteness, below we discuss only the region x < 0.5. At x > <5 
the probability of interaction of the qgq system with the medium (it is of the order of 
nas(pd, x)L) becomes small. Thus, it is clear that in the developed diffusion regime the 
spectrum is dominated by the N — 1 scattering. It is surprising that this turns out to be 
in apparent contradiction with prediction of the LLA. The LLA spectrum can be obtained 
using in (Bf) the oscillator Green's function. For zero parton masses it gives 



dP 

dx 



2G{x) 



7T 



l-J z 

Re J dz J d£ 



Q 2 



2G{x) 



cos 2 f2£ 



In I cosfiLl 



7T 



where G(x) = a s C F [l — x + x 2 /2]/x . This spectrum has been derived in 
\flL\ ~ 1 at x ~ 5. For the diffusion regime from ( |I4] ) one gets 



dP_ 

dx 



LLA 



z>5 



G(x)C 2 n 2 L 4 
8nE 2 x 2 (l - x) 



2 ' 



(14) 
Note that 

(15) 



Since the right-hand side of ([15]) oc n 2 it is clear that it corresponds to the N = 2 term. 
Thus one sees that the N = 1 contribution is simply absent in the LLA. 

The fact that the LLA fails in the diffusion regime can be directly seen from calculation 
of the N = 1 contribution. To obtain it one should use in ([]) the free Green's function 
(E). Then in the massless limit 



da* f H f (x,z) 



dx 



gives 
G(x)A(x)M(x 



N=l 



2tt 



Im 



d£ 



dp p C 2 (p) exp 



' iM[x)p 2 



(16) 



For C*2(p) = const the p 2 -integral in (^) has zero imaginary part, and the right-hand side 
of ( |16|) is also zero. On the other hand, using ([TT| ) one gets from flTB] ) 



da* f H f (x,z) 



dx 



N=l 



Then (|) yields 



dP_ 

dx 



N=l 



a 2 7iC T C F G(x)A(x)z 
AEx(l - x) 

a 2 7iC T C F G(x)A(x)nL 2 
8Ex(l-x) ' 



Let us see why the LLA fails in momentum representation in which (01) reads 



da* H f (x,z) _ a 2 C T C F A(x) 



Re / dpdq 



[tf *(p, x) - ^*(p - q, x)]* m (p, x, z) 



dx (2tt) 2 
In the massless limit from dl9|) one can obtain 

a 2 C T C F G{x)A{x) 



2\2 



d^ f H f (x,z) 



dx 



N=l 



2tt 



(q 2 + p 



d P 2 d q 2 F ^ q) 



(17) 



18) 



(19) 



2\2 



F(p, q) = Re — 

p2 



exp 



izp 



2M{x) 



2tt 



(q 2 + V 



q(q - pj 



p) 2 ' 



(20) 
(21) 
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where is the angle between q and p. The logarithmic situation with dominance of 
q 2 <C p 2 would correspond to F(p, q) oc q 2 at q 2 <C p 2 . However, the azimuthal integral 
in (PH) equals 2^9{q 2 — p 2 ), and the process is dominated by hard t-channel exchanges 
with q 2 > p 2 ~ 2M(x)/z. After integrating over p 2 and q 2 in (^) one reproduces (|17|). 

It must be emphasized that the LLA fails only in the diffusion regime. But it is a 
good approximation in the infinite medium regime when \l/ m falls off rapidly at the scale 
much smaller than Pd{L). It is also worth noting that the boundary (|T3|) beyond which 
the diffusion regime occurs is obtained assuming that in the infinite medium regime LPM 
suppression is strong (it means that p e ff(x ~ (5) < 1/e). It is possible that there the 
Bethe-Heitler situation takes place. One can easily show that in this case 5 ~ Lp 2 /2E. 
Thus, in general, the diffusion regime occurs for the gluons with energy u > co> cr , where 

f nC 3 L 2 Lp 2 \ 

Wcr ~ max — - — 1~2~) • ( 22 ) 
Let us now discuss the energy loss. It can be written as 

AE = [ duu— + / duuj— . (23) 

One can show that the first term in (^) does not depend on energy, and is of the order 
of AE BDMS (HD for both the LPM and Bethe-Heitler situations. At E — > oo the energy 
loss is dominated by the second term in (|23|) which grows logarithmically with E. Then, 
using flTBD to the logarithmic accuracy one can obtain in the high-energy limit 

. „ C F a s L 2 fi 2 E 

AE = ^—^\og — . (24) 

Here we have used a^7rCj?CTA(0)n/2/i 2 = 1/A S . Note that since L ^> from (|22|) it 
follows that always u cr 3> /x. The qualitative estimates (including the region uj < u cr ) 
show that the appearance of u cr in the logarithm in ( 124] ) instead of fi in (0) for RHIC 
conditions (L ~ 4 fm) can suppresses the energy loss at E ~ 10 GeV by a factor of 
~ 0.5. For SPS conditions (L ~ 2 fm) the suppression is not strong (~ 0.7 — 0.8 at E ~ 5 
GeV). The above estimates are obtained for the plasma temperature T = 250 MeV. Note 



that the absence of cu cr in the GLV prediction (y) is connected with the neglect in [13[ of 
the mass effects in evaluating the phase factor which controls the interference for gluon 
emission from different points of the quark trajectory. 

The above analysis is valid for the gluon emission from a fast gluon as well. In this 



case in ([24] ) Cf should be replaced by 2Ca (here the factor 2 comes from symmetry of 
the spectrum with respect to change x <-» (1 — x)). 
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